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The pD atoms are studied in models of various realistic, popular NN potentials. The small 
energy shifts and decay widths of the atoms, which stem from the short-ranged strong interactions 
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between the antiproton and deuteron, are evaluated in a well-established, accurate approach based 



on the Sturmian functions. The investigation reveals that none of the employed potentials, which 

X> ■ — 

^ , reproduce the NN scattering data quite well, is able to reproduce the experimental data of the 

energy shifts of the 2p pD atomic states. The energy shifts of the 2p pD atomic states are very 
sensitive to the NN strong interactions, hence the investigation of the pD atoms is expected to 
provide a good platform for refining the NN interaction, especially at zero energy. 
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The second simplest antiprotonic atom is the antiprotonic deuteron atom pD, consisting of an antiproton and a 

deuteron bound mainly by the Coulomb interaction but distorted by the short range strong interaction. The study of 

the pD atom is much later and less successful than for other exotic atoms like the protonium and pionium. Experiments 

in : were carried out at LEAR just in very recent years to study the properties of the pD atom Q,Q|- Even prior to the 
OO 

Q\ • experiments some theoretical works 



interactions. Recently, a theoretical work 



had been carried out to study the pD atomic states in simplified pD 
proposed a mechanism explaining the unexpected behavior, of the 
scattering lengths of NN and pD system, that the imaginary part of the scattering length does not increase with the 
size of the nucleus. 

In the theoretical sector, one needs to overcome at least two difficulties in the study of the pD atom. First, 
the interaction between the antiproton and the deuteron core should be derived from realistic NN interactions, for 
example, the Paris NN potentials 0,3,3, the Dover- Richard NN potentials I (DR1) and II (DR2) [ZZl , and the 
Kohno-Weise NN potential [12]. Even if a reliable pD interaction is in hands, the accurate evaluation of the energy 
shifts and decay widths (stemming for the strong pD interactions) and especially of the nuclear force distorted wave 
function of the atom is still a challenge. It should be pointed out that the methods employed in the works 
are not accurate enough for evaluating the wave functions of the pD atoms. 

In the present work we study the pD atom problem employing a properly adapted numerical method based on 

n, 

Sturmian functions 13]. The method accounts for both the strong short range nuclear potential (local and non-local) 
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and the long range Coulomb force and provides directly the wave function of the pD system with complex eigenvalues 
E = En — i^. The protonium and pionium problems have been successfully investigated [14I [15I] in the numerical 
approach. The numerical method is much more powerful, accurate and much easier to use than all other methods 
applied to the exotic atom problem in history. The pD interactions in the work are derived from various realistic NN 
potential, which is state-dependent. The work is organized as follows. The pD interactions are expressed in Sec. II 
in terms of the NN interactions. In Sec. Ill the energy shifts and decay widths of the Is and 2p pD atomic states 
are evaluated. Discussions and conclusions are given in Sec. Ill, too. 



II. pD INTERACTIONS IN TERMS OF NN POTENTIALS 



We start from the Schrodinger equation of the antiproton-deuteron system in coordinate space 
( P 2 P 2 \ 



2M„ 2M 



A 

where A and p are the Jacobi coordinates of the system, defined as 

A = r3 2 — ' P = r 2- r i (2) 

M p = M/2 and M\ — 2M/3 are the reduced masses. Here we have assigned, for simplicity, the proton and neutron 
the same mass M . Eq. |TJ can be expressed in the form, where the strong interaction is expressed in the isospin basis, 

{^A + VS + VC )^- E ^ (3) 
where V$ and Vq stand for the nuclear interaction and Coulomb force, respectively, and take the forms 



V S = V° N (r 2 fi) + \{V£ N (r 3 n) + V-£ N (r 3 - f 2 )} + - n) + V± N (r 3 - r 2 )] (4) 



Vc = \ [Vc {r 3 -n) + V c (r 3 - r 2 )] (5) 

V a and V 1 in eqs. ([!]) and (U|) are the isospin and 1 nuclear interactions, respectively. Note that we have assigned 
r*i2 as the relative coordinate of the deuteron core. 

One may express the interactions Vc and V$ in eq. ([U eq. ([5]) in terms of the interactions of certain NN states. 
In the I J MLS) basis of the pD states 

\JMLS) = \{L p ® L X ) L ® (S 12 ® S 3 ) S ) JM (6) 

we derive 

(H + W C (X, P ) + V° N ( P ) + W S (X, p)) *(A, P ) = E*{\, p) (7) 

with 

P 2 P 2 

H = — p — + — (8) 

2M p 2 Ma V ; 
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Wc and Ws in eq. ([7]) are respectively the Coulomb force and strong interaction between the antiproton and deuteron, 
and the interaction between the proton and neutron in the deuteron core. Wc and Ws are derived explicitly as 

(9) 



W c {\p) = \j dxV c (r 13 ) 



with 



W s (\,p) = \ [ dxJ2 (P\Q) (Q\ Vw N (ri 3 ) \Q') (Q'\P' 
J - 1 Q,Q' 



ri3 



\n~r 3 \ = (\ 2 + p 2 /4-Xp: 



1/2 



where x = cos 9 with 9 being the angle between A and p. In eq. (p~0|) \P) = \ JMLS) and \P') 
defined in eq. © while the states | Q) and | Q') are 

\Q) = \(L a <g> S 13 )j„ ® (L 7 <g> 5 2 )j 7 > JM 

\Q') = \(L' a ®S 13 )j„®(L 7 ®S 2 ) Jy ) JM 
Here <r and 7 are also the Jacobi coordinates of the system, defined as 



7 = r 2 



r 3 - n 



(10) 

(11) 

(12) 

\JML'S) are as 

(13) 
(14) 

(15) 



So defined the states | Q) and | Q') is based on the consideration that the NN interactions can be easily expressed in 
the I J a M a L a S\ 3 ) basis of the NN states. Note that (P\Q) depends on not only the quantum numbers of the states |P) 
and \Q), but also A, p and the angle 9 between A and p resulting from the projection of the orbital angular momenta 
between different Jacobi coordinates. We listed the integral kernels in eq. (fT0|) . J^q q> (P\Q) (Q\ ^(^13) \Q') {Q'\P')i 
for the lowest pD states in the approximation that the deuteron core is assumed in the S-state, as follows: 
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P') = 
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(16) 
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where \P) = \JMLS) and \P') = \ JML'S) are the pD atomic states. Both the pD and NN states in eq. (QU) 
are labelled as 2S+1 Lj with S, L and J being respectively the total spin, total orbital angular momentum and total 
angular momentum. The potentials Vj^ N , being functions of r\3 = \/A 2 + p 2 /4 — pXx, stand for the NN interactions 
for various NN states as indicated in the brackets. 

The F\, F2 and F$ in eq. (fT6|) are functions of only A and p, taking the forms 



P<2\ 



4A 
3(0' 



p > 2A 



(17) 



F 2 = 



1 - ) , p < 2A 



0. 



p > 2A 



(18) 



F 3 



9j plfl _3 3 J-\ 
2 J 2 ' 2 ' 4A 2 '' 
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V p 
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8A 



3p J 
64A :i 



p < 2A 
p > 2A 



(19) 



where 2-F ll (a, /3, 7, x) is the hypergeometric function and Artanh(x) the inverses hyperbolic tangent function. 



III. ENERGY SHIFTS AND DECAY WIDTHS OF pD ATOMS 

It is not a simple problem to accurately evaluate the energy shifts and decay widths, especially wave functions of 
exotic atoms like protonium, pionium and antiproton-deuteron atoms, which are mainly bound by the Coulomb force, 
but also effected by the short range strong interaction. In this work we s tudy the pD atoms in the Sturmian function 
approach which has been successfully applied to our previous works Q, (3- Employed for the NN interactions are 
various realistic NN potentials, namely, the Paris NN potentials of the 1994 version (Paris84), 1998 version (Paris98) 
and 2004 version (Paris04), the Dover- Richard NN potentials I (DR1) and II (DR2), and the Kohno-Weise NN 
potential (KW). In this preliminary work, we just limit our study to the approximation of undistorted deuteron core. 
However, one may see that the main conclusions of the work are free of this approximation. 

Shown in Table I are the energy shifts and decay widths, which stem from the Paris98, DR2 and KW NN interac- 
tions, in the approximation of undistorted deuteron core. The theoretical results for other interactions like Paris84, 
Paris04 and DR1 are quite similar to the ones listed in Table I. The wave function of the undistorted deuteron core is 
evaluated in the Bonn OBEPQ potential 16]. It is found that the theoretical results for the Is pD atomic states are 
more or less the same by all the employed NN potentials. The predicted energy shifts are roughly as twice large as 
the experimental data. However, one may expect that the predictions of the potentials in question could be improved 
to some extent by solving the pD dynamical equation in eq. ([7]) without any approximation. A better treatment of 
the deuteron core will yield lower Is pD atomic states, hence smaller energy shifts. The theoretical results for the 
decay widths of the Is pD atoms are also larger than the experimental data though not as far from the data as for 
the energy shifts. The predictions for the decay widths are also expected to be improved by treating the deuteron 
core more properly. 

The theoretical predictions for the energy shifts of the 2p pD atomic states are totally out of line for all the NN 
potentials employed. The experimental data show that the averaged energy level of the 2p pD atoms is pushed up 
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-2445 
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-2673 
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2450 






4 SD 3 / 2 


-2680 
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-2668 
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-2503 


2469 






2 p 
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99 


657 
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265 


402 


47 


846 


101 


785 






2 -P 3 /2 


-128 


515 


14 


897 


98 


643 






4 PF 3/2 


282 


477 


21 


887 


97 


648 






4 PF 5/2 


244 


814 


21 


877 


101 


660 






AEi s , Tig 


-2602 


2475 


-2670 


2387 


-2494 


2463 


-1050 ± 250[1J 


1100 ± 750[1] 
2270 ± 260[2] 


AE2 P , T 2p 


124 


602 


22 


883 


99 


668 


-243 ± 26[2] 


489 ± 30[2] 



TABLE I: The energy shifts AE and decay widths of the Is and 2p antiproton-deuteron atomic states in the approximation 
of undistorted deuteron core. The minus sign of the energy shifts means that the strong interaction is repulsive. The units are 
eV and meV for Is and 2p states, respectively. 

by the strong interaction, the same as for the Is pD atoms, but the theoretical results uniquely show the averaged 
energy level shifting down. It is unlikely to improve, by treating the deuteron core more accurately, the theoretical 
predictions of the NN potentials in question for the 2p pD energy shifts since a more accurate treatment of the 
deuteron core will lead to deeper 2p pD atomic states. 

All the NN potentials employed in the work reproduce AW scattering data reasonably, but badly fail to reproduce 
the energy shifts of the 2p pD atoms. The investigation of the pD atoms may provide a good platform for refining 
the AW interaction, especially at zero energy since the energy shifts of the 2p pD atomic states are very sensitive to 
the A^ A~ strong interactions. 

The research here is just a preliminary work, where a frozen, S-state deuteron is employed. The work may be 
improved at two steps, considering that the numerical evaluation is time-consuming. One may, at the first step, solve 
the pD dynamical equation in eq. ([7]) by expanding the pD wave function in a bi-wave basis of the Sturmian functions, 
where a realistic nucleon-nucleon potential is employed but the deuteron core is assumed to be at the S-state. Such 
an evaluation is still manageable at a personnel computer but it may take a week or longer. We may compare the 
results of the improved work with the results here to figure out how important an unfrozen deuteron core is. 

One may also consider, at the second step, to solve the pD dynamical equation in eq. ([7|) by expanding the pD wave 
function in a bi-wave basis of the Sturmian functions without any approximation, where realistic nucleon-nucleon and 
nucleon-antinucleon potentials are employed and the deuteron core is allowed to be at both the S- and D- waves. It is 
certain that the numerical calculation will take longer time but, anyway, we will do it after we complete the first-step 
improvement. 
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